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Abstract: We show how quark-disconnected and quark-connected contributions to 
hadronic n-point functions can be written as independent correlators for which one can 
derive expressions in partially quenched chiral effective theory. As an example we apply 
the idea to the case of the hadronic vacuum polarisation. In particular, we consider 
the cases of the Nf = 2 theory without and with a partially quenched strange quark 
and also the Nf = 2 + 1 theory. In the latter two cases a parameter-free prediction 
for the disconnected contribution at NLO in the effective theory is given. Finally we 
show how twisted boundary conditions can then be used in lattice QCD to improve 
the q 2 resolution in the connected contributions even when flavour singlet operators are 
considered. 
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1 Introduction 



Hadronic correlation functions which involve quark disconnected contributions are no- 
toriously hard to compute in lattice QCD [U[2]. Such contributions however appear 
in several quantities related to properties of flavour singlet particles (e.g. the 77 and 
7/ mesons) or to matrix elements of flavour singlet operators (as for the strangeness 
content of the nucleon), or to electro-magnetic interactions (e.g. the hadronic contribu- 
tion to the anomalous magnetic moment of the muon or to the nucleon electric dipole 
moment). There is no conceptual difficulty in treating disconnected quark diagrams but 
the computational effort is immense compared to quark-connected contributions. The 
disconnected part is therefore often neglected, not always providing solid arguments 
that the systematic uncertainty introduced in this way is under control. 

Here we present a method that allows to predict the quark disconnected contribu- 
tion to correlation functions in (partially quenched) chiral perturbation theory [3HZ]- 
We first give the general argument and then present an explicit example by deriving 
predictions for the magnitude of the quark disconnected contribution to the hadronic 
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vacuum polarisation of the photon for Nf = 2 dynamical flavours, for Nf = 2 dynam- 
ical flavours with a quenched strange quark and also for the case of Nf = 2 + 1 active 
flavours. The vacuum polarisation is the main ingredient in computations of the lead- 
ing hadronic contribution to the anomalous magnetic moment of the muon and hence 
of relevance to precision tests of the Standard Model (see Ref. [8] for a review of the 
subject). 

We will stay in the Euclidean continuum and infinite volume for most of the discus- 
sion. While having applications to lattice QCD in mind, reference to this regularisation 
is only made where we think that it helps the better understanding of our arguments. 
We will discuss however the finite volume case in order to show how this technique 
allows for using partial twisting [91-TH] also for form-factors and polarisations involving 
flavour-diagonal operators. 

2 General argument 

We consider a n-point fermionic correlation function in QCD with Nf dynamical flavours 
(not necessarily degenerate). In general several Wick contractions contribute to the 
correlator. By introducing valence quarks which are degenerate with the dynamical 
flavours, each Wick contraction can be rewritten in terms of a single fermionic cor- 
relation function defined in an un-physical theory. The physical result is recovered 
by summing over the correlation functions in the un-physical, partially quenched, the- 
ory [5H7]. In particular, the approach can be used to separate the contributions from 
quark disconnected diagrams from those coming from quark connected diagrams. These 
contributions taken on their own are un-physical, it is therefore natural that in order to 
define them as correlation functions one has to resort to un-physical theories. Investi- 
gations of quark-disconnected diagrams along these lines can also be found in Ref. [7]. 

The number of valence quarks N v (i) degenerate with the ith dynamical flavour, 
which has to be introduced depends on the particular correlation function. Given a 
n-point fermionic correlation function, N v (i) is related in an obvious way to the largest 
number N™ x (i) of disconnected quark loops involving the ith flavour, which can appear 
as all possible Wick contractions are considered. In particular N v (i) = Nj^fy) — 1 at 
most, in such a way that a different flavour can be attached to each disconnected quark 
loop. In this approach the partially quenched theory, in which each Wick contraction 
of the original correlator can be written as an independent correlation function, is not 
unique but rather depends on the specific n-point function we started from. 

Partially quenched chiral perturbation theory (PQxPT) [3HZ], which is an extension 
of chiral perturbation theory [31 0], provides an asymptotic low energy description of 
partially quenched QCD (PQQCD) and can therefore be used to obtain predictions 
for, and algebraic relations amongst, the disconnected and the connected part of a 
fermionic correlation function. There is a rather large body of literature on PQ^PT 
and PQQCD, we found the reviews in Ref. [T5] and Ref. [16] very clear and useful. 
In short, following Ref. [15], PQQCD can be formulated in terms of a local theory by 
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introducing a commuting spin-1/2 field, a ghost [17], labelled by q for each valence 
quark q and by extending the fermionic QCD action by including terms of the form 
q{p+m q )q-\-q^ {fi + m q )q, which violate the spin statistics theorem. As long as all masses 
are positive, the determinant produced by the integral over the valence quark fields is 
cancelled by the corresponding integral over the ghost fields and the QCD partition 
function is reproduced. The advantage is that this formulation provides field-theoretic 
expressions for partially quenched correlation functions. As discussed in Refs. |15|I16|. 
a low energy description of the theory (PQxPT) is obtained by constructing a chiral 
Lagrangian encoding the apparent chiral symmetry group of the extended action, which 
is a SU(N } + ES N v (i) | ES JV«(«)) graded group. 

For the application discussed here the strategy can be exemplified for the case of a 
meson 2-point function constructed of flavour-diagonal quark-bilinears, 

Cqcd(v,x) = (9(y)r'g(y)g(x)rg(x)) QCD 

q=u,d,s,... 

(2.1) 

= - ]T (Tr{S q (x,y)T'S q (y,x)T}) QCD 

q=u,d,s,... 
q=u,d,s,... 

where S q is the propagator of the quark field q and rW may contain Dirac- as well 
as colour-structures. The trace is over Dirac- and colour-indices. We will rewrite the 
correlation function Cqcd(u,x) as the sum of two new correlation functions: one being 
equivalent to the quark-connected contribution and the other being equivalent to the 
quark-disconnected contribution. To this end for each quark q we add a valence quark 
q v which is mass degenerate to it, together with the corresponding ghost field. The 
correlation Cqcd(u,x) can be rewritten as 

Cqcd{v,x) = ^2 (q(y)r'q v (y)q v (x)rq(x)) pQQCD (2.2) 

q=u,d,s,... 

+ Yl (q(y)r'q(y)q' v (x)Tq v ^ 



'PQQCD 



q=u,d,s, 
iConn /„, \ i /-(Disc 



= CpQQCD(y> x) + Cpqqcd(2/; x) . 

The Wick contractions of the first correlator on the r.h.s. of the first equation lead to 
quark connected diagrams only, while those for the second correlator produce quark- 
disconnected diagrams only (see illustrations in figure [1]) . If we imagine for a moment 
having regularised the theory on the lattice it is clear that the above equality holds 
non- perturb at ively, since on each gauge configuration the quark propagators and the 
fermionic determinant (and therefore the weight of the configuration in the path integral) 
are the same in the two theories. In the following we will show how ^pqqcd^' x ) anc ^ 
CpQQ CI) (y,x) can be expressed in a suitable chiral effective theory, depending on the 
flavours entering the sum in Eq. 12.11 
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(a) (b) 
Figure l: Wick contractions: (a) connected and (b) disconnected diagram. 

3 Hadronic vacuum polarisation 

The Euclidean hadronic vacuum polarisation (VP) tensor is defined as 

ng^te) = i J dW^J^xJJ^O)) = iC^(q) , (3.1) 

/AT \ J 

where J M / (x) = Yl QqQ( x )l^Q{ x )- For Nf = 2, q = (u,d) and Q q = (2/3, —1/3) and 

9=1 

for Nf = 2 + 1, q = (u,d,s) and Q q = (2/3,-1/3,-1/3). Euclidean invariance and 
current conservation imply 

ng^Gz) = (q,qu ~ 9^q 2 )rt N f\q 2 ) . (3.2) 

For space-like momenta, the relation between Ilgr (g 2 ) and the lowest order hadronic 
contribution a^ LO to the anomalous magnetic moment of the muon has been derived in 
Ref. [T8j[19] and reads (suppressing the index Nf) 

^ 2 

with 



C° = ("YI tf/foW). (3-3) 



, m 2 q 2 Z 3 (l-q 2 Z) q 2 - Jq A + 4m 2 q 2 

f^l) = 1 _|_ a 272 , ^ = \ 2 2 • ( 3 - 4 ) 

and I% 2 ) = 4tt 2 [n(g 2 ) - 11(0)]. We consider the iso-scalar meson two-point function 
C iiu {q) in Eq. (|3.ip . The Wick contractions for that correlator lead to quark diagrams 
of type (a) and (b) as illustrated in figure [TJ By considering in some detail the two 
flavour case we will first separate the disconnected contributions from the connected 
ones in the way discussed in the previous section and then set-up a computation in the 
resulting PQxPT framework. 

According to the discussion above we would need two valence quarks, one degener- 
ate with the u quark and one degenerate with the d quark. However, as we will always 
assume iso-spin to be an exact symmetry, it is enough to introduce one valence quark, 
which we call r, with the corresponding ghost r g . We therefore have an SU(3\1) chiral 
group. The correlation in Eq. (|3.1j) can then be decomposed as 



= /^xe^^(j7(x)j«(0)) + i(j*(^(0)) + 

(3.5) 
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where jfi iq2 (x) = qi(x) / y^ l q 2 (x). The first two correlators on the r.h.s. are connected, 
whereas the last three represent the disconnected contributions to the hadronic VP 
tensor. It is convenient to cast the quark-fields into a four-component vector ^ T = 
(u, d, r, r g ) and introduce the generators T a , a = 1, . . . , 15 of the graded group SU(3\1). 
We use the conventions also employed in [20] , 

(3.6) 



T a = (T a ) t , Str{T a }=0, Str{T a T 6 } = ^g ab , o = l,...,15, 
with the super-trace Str {^4} = An + A 22 + ^33 — A44. The matrix g ab reads 



fl 



\ 



-a 2 



-a 2 



-a 2 



V 



1 



9-14 



(3.7) 



•v 



} 15 



where a 2 is the second Pauli matrix. T 1 , . . . , T 8 are the generators of the SU (3) sub- 
group that acts on the sea and valence components, T 9 , . . . ,T 14 mix the quark with 
the ghost components, T 15 is the diagonal matrix diag(l, 1, l,3)/(2v 3). We also add 
T° with Str{T } = l/y/2 which is proportional to the unit matrix in order to describe 
iso-scalar interactions. For our choice of the generators T°, . . . , T 15 , the currents in 
Eq. (|3.6p can then be rewritten as 

= ^{x)^(T 4 + iT 5 )^(x) : 

= ^(x) 7m (T 6 + zT 7 )^(x), 
: ij(x) lfl (T 6 -iT 7 )ij(x), 



3j W 

f;{x) 



j™(x) 



ij(x) lfl \y2T° - ^T 15 + ±T 8 + T 3 J i;(x 
${x) llx (y2T° - ^T 15 + ^T 8 - T 3 ) i{j(x 
ij(x) 7 » (V2T 



(3. 



1 7^15 

V3 



^)^(X) 



This form is more suited for the PQ%PT computation, as it will become clear in the 
following. 

Before concluding this section we note that by using again iso-spin symmetry we 
could have separated the disconnected and connected parts in the correlator above with- 
out introducing any additional valence quark. We have used this completely equivalent 
approach in Ref. [21]. Here however, we preferred to introduce a graded flavour group 
to provide an example for how one has to proceed for the the more general but also 
more complicated case of 2 + 1 flavours. 
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4 PQxPT for the connected and disconnected parts of the hadronic 
VP 



In this section we briefly introduce those parts of the (partially quenched) chiral La- 
grangian up to 0(p A ) which contribute to the connected as well as to the disconnected 
piece of the VP [3 , 4J . We took care that the discussion applies to any choice for the 
(graded) symmetry group. In the next section we will then present results for the cases 
of SU(3\1), SU(4:\2), and SU(4:\1) flavour groups. These are the relevant symmetry 
groups for the description of the contributions to the VP in the Nf = 2-theory without 
and with a quenched strange quark and for the Nf = 2 + 1-theory, respectively. 



4.1 0(p 2 )- Lagrangian 

For a generic graded flavour group, the leading order chiral Lagrangian is [3H7] 

£( 2 ) = ^Str {DyUDptf} - \BF 2 $ti {mU ] + M ] u} , (4.1) 

where D^U = d^U + iv — iU v ^ provides the coupling of the meson field to an external 
vector source Vn and where 

M = diag (m x ,...,m N} , mf cacc , . . . , m^ cacc , mf host , . . . , m^ host ) , (4.2) 

contains the dynamical, valence-, and ghost-quark masses. By chirally expanding 
U = exp \2i &-jjr ^, where T a are the N = (Nf + 2N) 2 — 1 generators of the corre- 
sponding SU(Nf + N\N) graded symmetry group (with the generators normalised as 
in the previous section) and where the 4> a are the Goldstone-boson/fermion fields, we 
arrive at the 0(p 2 ) -expression 

C {2) = \ g ab d^ a d^ b + 2£M ffl W + 42 , (4.3) 

where M ab = Str {MT a T b j. With = u°T a , the part of the interaction Lagrangian 
that is of relevance here is 

N N N 

42 = - E c k bc 9 a %rv b ^ c +^EE c^c^v^^i , m 

k=l k=X 1=1 



a) b) 

where 

TV 

C a bc = -2i Str { [T b , T c ] T fc } g ka , (4.5) 

k=l 

are the structure constants of the underlying (graded) symmetry group with the 
(anti-) commutator 

[T a ,T b ] = T a T b — (—) Va Vb T b T a . (4-6) 
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(a) (b) (c) 

Figure 2: Contributing vertices in the effective theory. 







x 1 


x 2 


SU(2), 


SU(3\1), SU(A\2) 


-4/i 2 


h s 


SU(3), 


SU(4 1) 


Lio + 2Hi 





Table l: Gasser-Leutwyler coefficients. 



The r\ a (and correspondingly for %) are 1 if T a mixes valence or sea quarks with a ghost 
and otherwise (cf. [22] ). The Feynman-rules for the vertices illustrated in figure [2] (a) 
and (b) are then determined as 



' ba „/ 

"a* 



a) -(C^pv + C* 

n . ' (4.7) 



k,l=l 

with c a 6c = c^V*- 

4.2 0(p 4 )-Lagrangian 

The relevant terms of the 0(p 4 ) Lagrangian [3j[4"ll231[24"]. i.e. those parts that have a 
non-zero matrix element between single external vector-sources, are 

£W = ZiStr{VV} + ^Str{^}Str{^}, (4.8) 

with 

VfW = Q^Vy - O^VfJ, , (4.9) 

and where is the trace-less part — ^- Strjv^j,}. As summarised in tabled] the 
coefficient X\ is the shorthand notation for the Gasser-Leutwyler low energy constants 
of the underlying symmetry group (see for example [3jSl[25] ) and X 2 parameterises 
the effective dynamics of flavour-diagonal contributions. The vertices corresponding to 
these Lagrangians have the following form (cf. figure [2] (c)): 



c) i 2 Xi g ab (v^p 2 - PfjPu) for a, b > , 
c) i 2 X 2 (rjfiuP 2 ~ PfiPu) for a = b = . 



(4.10) 
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Figure 3: Feynman diagrams contributing to the hadronic VP at 1-loop level: (a) 
unitary, (b) tadpole, (c) 0(p 4 )-insertion 

4.3 Contributing diagrams 

At NLO only the diagrams in figure [3] (a) and [3] (b) contribute dynamically to the VP 
and the diagram in figure (c) provides counter terms. Let tp be a flavour vector of 
Nf dynamical quarks, N partially quenched valence quarks and -/V corresponding ghost- 
quarks. The expression in the effective theory for the Fourier transform of the correlator 
i('ip(x)^f ll T a ilj(x) ip(0)^f l/ T b ip(0)} at NLO in chiral perturbation theory is of the form 

rr>/) ( q , /x) = n (^ unit - ( q ) + nfr b ^- ( q ) + u^ co ^- ( M ) , (4.11) 

where \x is the renormalisation scale. The first term on the r.h.s. corresponds to the 
unitary diagram, 

U M,nnit. (q) = ^C. aj C k bl (4.12) 

x / 0p{q + Zk),{ q + 2k) u & k (k 2 )Gi l (( q -kf), 

where the overall factor | is a symmetry factor. For the second term, the tadpole- 
diagram (b), one derives 

n£j>)tad P . = in^cfcp J '-0p(P k (k 2 ). (4.13) 

For the SU(3\1) theory the meson propag ators G ab (k 2 ) can be found e.g. in Ref. |2U| 
and we provide the corresponding expressions for SU(4\1) and SU(4\2) in appendix [Al 
The last term, n^^' count ' (;u), contains the counter terms. We compute U^ b \ q , fi) in 
dimensional regularisation as explained in detail in Ref. |26j where also the solutions for 
the loop integrals have been derived and reduced to the integral £?2i(<? 2 , M 2 ) defined in 
appendix [Bj 

5 Results and applications to Lattice QCD 

Here we present the results for the VP at NLO in chiral perturbation theory for various 
underlying graded flavour-symmetry groups. In particular, we consider the theories 
with Nf = 2 flavours without or with an additional quenched strange quark and also 
the Nf = 2 + 1 theory. In each case we present the expressions in the effective theory for 
the full VP and also for the contributions from quark-connected and quark-disconnected 
diagrams. 
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5.1 N f = 2 in SU{3\1) PQxPT 



For the theory with two dynamical and degenerate light quarks only we found the 
following expressions for the VP: 

nSV) =- (^^)+lh s + iAB 21 ^\ q \Mt)), 

ng^ 2 ) = - ¥(a<»W00 + i4B 2l (^, q 2 ,M 2 )) , (5.1) 

ngliV) = i (a (3|1) M-2^ +^ 21 (^W,m2) 



where 



A( 3 l 1 )(//) = -8/ l2 (^). 



(5.2) 



The integral E>2i(fJ. 2 , q 2 , M 2 ) is defined in appendix [Bj All contributions are param- 
eterised in terms of the low-energy constant h% and the pion mass M T . The full ex- 
pression as well as the disconnected piece also depend on the parameter h s . This is a 
peculiarity of the two-flavour theory and we will see in the next sub-sections that the 
corresponding expressions in the presence of a strange valence quark do not depend on 
this additional parameter. Hence, a prediction of the quark-disconnected diagram is 
not possible without the knowledge of h 2 and h s . These parameters however, can in 
principle be determined from simulations of lattice QCD and h s can also be obtained 
by matching to the expressions in the next section. 

The application we have in mind when discussing the VP is the leading hadronic 
contribution to the muon anomalous magnetic moment defined in Eq. (j3.3l) . As discussed 
in section[3l a^ LO depends on Tl(q 2 ) = 47r 2 [n(q 2 ) — 11(0)] . We observe that all reference 
to the low-energy constants disappears in this difference at NLO in the effective theory 
and we find |21] 

npisc(g 2 ) _ i . . 

At this order in the effective theory it is therefore sufficient to compute the quark- 
connected piece contributing to and then correct for the quark-disconnected piece 
using the above relation which predicts a 10% negative shift for all values of the momen- 
tum (note in this context the discussion in section I5.4|) . The disconnected contribution 
at NLO turns out to have the same momentum and quark-mass dependence as the 
connected piece. 



5.2 Nf = 2 plus a quenched strange quark in SU(A\2) PQxPT 

For the theory with two dynamical and degenerate light quarks and one quenched 
strange quark one needs to consider SU(4\2) PQxPT, where we found the following 
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expressions for the vacuum polarisation: 
nffiV) =A( 4 l 2 )( / u)-4i( B 21 ^ 2 ,q 2 ,M 2 ) - \B 2 ^\ q \Ml) + ±B 21 (v\q\Ml)) , 

nS 2 n ) n(« 2 )=A (4|2) (^)-4z( f B 21 ^,q\Ml) + \B 2X {p? , q 2 , M*j) , 

nKiV) = -4i(-iS 21 ( M 2 , g 2 ,M|) - \B 21 ^\q\Ml) + lB 2l {n\ q \Ml)) , 

(5.4) 

where 

A( 4 l 2 )(^) = 8/ l2 (/i). (5.5) 

The full expression and the connected piece depend on the low energy constant h 2 
and the pion and kaon mass. An additional dependence on the mass of the strange- 
quark enters in terms of the mass M 2 S = 2Bm s of an un-physical meson made by 
the two quenched strange quarks. It is well known that in the Nf = 2 + 1-theory 
the quark-disconnected contribution to the hadronic VP vanishes in the limit of equal 
quark-masses. This SU(3) symmetry prohibits the presence of low-energy constants in 
the expression for the quark-disconnected diagrams at NLO and this is indeed what we 
found when computing H^l^( q 2 ). At NLO in the effective theory we therefore provide 
an entirely parameter-free prediction of the quark-disconnected diagram. Moreover, 
as in the case of the Nf = 2-theory, also II ( q 2 ) is free of low-energy constants and a 
parameter-free prediction for the ratio of the quark-disconnected contribution to the 
quark-connected contribution can be made. As can be seen in figured] (a), in the limits 
Mk —> oo and Mk = M n , the above formulae reproduce the results for the Nf = 2 
theory and the vanishing of the disconnected piece in the Mk = M^-limit, respectively. 

5.3 N f = 2 + 1 in SU(4\1) PQxPT 

For the theory with two dynamical and degenerate light quarks and one dynamical 
strange quark we found the following expressions for the vacuum polarisation: 

IlffiV) =A( 4 ID( / u)-4i( B 21 ^ 2 ,q 2 ,M 2 ) + B 21 (» 2 ,q 2 ,M 2 K ) ) , 

n t :t ) „(9 2 )=A (4|1) ^)-4z( f B 21 ^,q\Ml) + \B 21 ^ 2 , q 2 , M 2 SS ) + |-B 2 i(/i 2 , q 2 , M K )j , 

4*(-i J B 21 ( M W,M 2 ) - lB 21 ^ 2 ,q 2 ,M 2 s ) + §5 21 (/x 2 , q 2 , M 2 K )) , 

(5.6) 

A( 4 l 1 ) = -2(L 10 (^) + 2if 1 ( / ,)) . (5J) 

As in the previous sub-section, the full expression and the connected piece depend on 
low-energy constants, this time Lio and Hi, and the pion and kaon mass. Also here a 
dependence on the mass of the strange-quark enters in terms of the mass M 2 S = 2Bm s of 



O 2 ) 



where 
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Figure 4: The plots illustrate the ratio of the contributions from quark-disconnected 
and quark-connected diagrams to the VP as a function of the momentum q 2 at NLO in 
chiral perturbation theory for the cases (a) Nj = 2 with a quenched strange quark, (b) 
Nf = 2 + 1. In both cases the lines are from top to bottom: the limit M n = Mk, then 
fixed M K = 495MeV and M w = 400, 300, 200, 139MeV. The bottom most-line at -1/10 
is the result for the Nf = 2 theory. 
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an un-physical meson made of a dynamical and a quenched strange quark. In the same 
way as in the previous section SU(3) symmetry prohibits the presence of low-energy 
constants in the expression for the quark-disconnected diagram at NLO. At this order in 
the effective theory we are therefore able to provide an entirely parameter-free prediction 
of the quark-disconnected diagram. Again, tl(q 2 ) is free of low-energy constants and 
a parameter-free prediction for the ratio of the quark-disconnected contribution to the 
quark-connected contribution can be made. We illustrate this result in figure H](b). In 
the same way as for the Nf = 2 theory with a quenched strange quark, in the limits 
limits Mk —> oo and Mk = M n , the above formulae reproduce the results for the Nf = 2 
theory and the vanishing of the disconnected piece in the 5i7(3)-limit, respectively. As 
a further check, our formulae reproduce the result in Ref. [26J. From the comparison of 
the plots (a) and (b) we see little influence of the dynamical strange quark. 

5.4 Discussion 

Before we suggest applications to simulations of lattice QCD we would like to point out 
obvious limitations of the above formulae. While the underlying idea of the approach, 
i.e. Eq. (|2.2p . is exact in QCD and in the chiral effective theory, the results presented 
here are only valid up to NLO chiral perturbation theory. If at all, this approxima- 
tion only holds for small quark-masses and only for small values of the momentum. 
There are in fact doubts that the physical strange quark mass can be described reliably 
within chiral perturbation theory [27] . One might hope that cancellations in the ratio 
nDisc( ( ? 2 )/IIconn(<? 2 ) will improve the convergence of the expansion in the meson masses 
and the momentum. However, only expensive numerical simulations in lattice QCD and 
higher order computations in the effective theory will be able to tell. 

At least as important is the fact that instead of coupling to two pions, the photon 
will also couple to vector resonances. It is a priori not clear if at all or in which 
kinematical regime chiral perturbation theory parameterises these dynamics correctly 
in terms of the low-energy constants of the chiral effective Lagrangian. Aubin and 
Blum |28] for example found indications for vector-dominance in the analysis of their 
lattice QCD data for the hadronic VP. We expect sizeable corrections to the expressions 
found here. There do exist models, generally known as "resonace xPT " models, where 
the vector resonances are dynamical degrees of freedom (see e.g. |29j). Conceptually 
these models are not as solid as xPT, to the extent that no systematic power counting 
can be formulated. While clearly desirable, a study of how the vector degrees of freedom 
will modify the expressions is beyond the scope of the current work but in principle 
technically straight forward within the framework of resonace xPT, which is however 
subject to the theoretical issues mentioned above. 

Given these remarks, the results presented here can only provide an order-of- 
magnitude estimate and should be used with great care. In the absence of any other 
quantitative knowledge of the quark-disconnected diagrams, the results can be used in 
order to correct the lattice data for the connected contributions. 
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5.5 Applications to lattice QCD 

Partially twisted boundary conditions for the VP: Partially twisted boundary 
conditions [9l-ll4j by now have become a standard tool in lattice hadron phenomenology. 
For the boundary condition qi(xk + L) = e i " k ' i qi(x) for the valence quark- flavours q± 
and q2 of a pseudo-scalar meson of mass m in a finite lattice of spatial extent L, the 
dispersion relation takes the form E(0i,6 2 ) = sjm 2 + (Jh/L - 6> 2 /L) 2 0H2HH!, which 
shows how the periodicity of the fermionic fields can be modified in order to induce 
spacial momentum to hadrons. Besides an exponentially suppressed and computable 
finite size effect the choice of quark-field boundary conditions for processes with only 
one initial and/or final hadron state introduces no further systematics [Sj. The effect 
of partial twisting cancels however in flavour-neutral mesons, as is clear from the form 
of the above dispersion relation. Due to the flavour-diagonal structure of the electro- 
magnetic current, partial twisting would therefore naively have no effect on the hadronic 
VP. The prescription introduced here however, allows for isolating a flavour-off-diagonal 
connected contribution and thus for inducing arbitrary values of the momentum for the 
VP. In other words we can assign different twisting angles to the quarks q and q v in 
Eq. 12.21 which has the net effect of inducing momentum in the connected part. The 
effectiveness of this method in the numerical computation of the connected part of the 
hadronic VP will be illustrated in a forthcoming publication [30J. The relations in the 
chiral effective theory between the quark-connected and the quark-disconnected piece 
then allow to predict also the momentum-dependence of the disconnected piece. A 
very similar argument was already used in Refs. |31| l32j in order to justify the use of 
partially twisted boundary conditions for the pion's vector form factor. The exponential 
suppression in the volume of the isospin breaking introduced by partially twisting only 
one of the light quarks is also discussed in Refs. [31|132|. 

6 Conclusions and outlook 

In this paper we show how quark-disconnected contributions to hadronic n-point func- 
tions can be estimated within the framework of partially quenched chiral perturbation 
theory. As an example we derive predictions for the quark-disconnected contribution 
to the photon's vacuum polarisation for QCD with two degenerate light quarks, with 
and without a quenched strange quark, and also for the case of the three-flavour theory. 
In the presence of a quenched strange quark the prediction for the quark-disconnected 
contribution turns out to be parameter- free. 

The vacuum polarisation is currently computed by various lattice collaborations 
with the aim to make first-principles predictions for the leading hadronic contribu- 
tion to the muon anomalous moment [HI [281 [33] . Since precise computation of quark- 
disconnected correlators from first principles still remain numerically extremely expen- 
sive they are often neglected in such computations. The approach presented here will 
allow to improve estimates of the systematic error introduced in this way. 
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The method can be applied to other physical observables of interest that receive 
contributions from quark-disconnected diagrams, like the strange-quark content in the 
nucleon, the nucleon electric dipole moment or the pion scalar form-factor. The discus- 
sion presented here is limited to NLO chiral perturbation theory. As a next step we 
plan to include the vector degrees of freedom and it would also be important to extend 
the computation for the case of the vacuum polarisation to higher orders in the chiral 
expansion in order to assess how robust the predictions reported here are. The inclusion 
of the leading lattice artifacts in the effective theory is also straight forward within the 
framework set up by Sharpe and Singleton in Ref. |34j . 

A spin-off of the ideas presented here is a method that allows for projecting (the 
connected part of) correlators containing flavour-diagonal currents onto any desired mo- 
mentum using partially twisted boundary conditions. 

Acknowledgments: We thank Dalibor Djukanovic, Silvia Necco, Hartmut Wittig, and 
Steve Sharpe for illuminating discussions. 
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Appendices 



A Propagators in chiral perturbation theory 

Here we provide the definition of propagators in PQxPT with the graded symmetry 
groups SU(4\1) and S77(4|2). For simplicity all results are quoted in the limit where 
the Goldstone bosons and the corresponding ghosts are mass-degenerate. This limit is 
sufficient for the discussion in the main body of this paper. For completeness we note 
that the case SU (3\1) has been worked out nicely in [20] and for most of the discussion 
we adopt the conventions laid out there. 



A.l Propagators in SU(4\1) chiral perturbation theory 

We consider the case of the effective theory with Nf = 2 + 1 dynamical flavours with 
two degenerate light quarks of mass m, one strange quark of mass m s and one quenched 
r quark with m s = m r = m Tg . The dynamical (not necessarily physical) degrees of 
freedom in the effective theory have masses 

M 2 = 2Bm, M 2 K = B(m + m s ), M 2 ss = 2Bm s . (A.l) 

The propagator G l i{k 2 ) is defined as follows: The Lagrangian in Eq. (|4H) with JCQ = 
can be rewritten as 

£® = I <f* [d^ a 8^ b + M a W) + \h ab {Ml - M 2 K )r^ b , (A.2) 

where 

Ml for a = 1,2, 3, 



M 2 = { ^{Ml + M 2 K ) for a = 4, 5, 6, 7, 9, 10, 11, 12, 16, 17, 18, 19, (A.3) 
M 2 S for a = 8,13,14,15,20,21,22,23,24, 



and where for our choice of the 5 , ?7(4|l)-generators 

h ab = (^g a8 + g al5 -g a24 )(V2g b8 + g bl5 -g b24 ) 

l*b = _^ g aS + g al5 +g a24^^2 g bS +g bl5 +g b2A- ) 



Gi(k\M 2 ) = IF ^ y 
G 2 (k 2 ,M 2 ,M 2 ) - 



(A.4) 



In terms of these matrices the propagator is given by 

G ab (k 2 ) = ^G 1 (fc^M a 2 ) + I/ afe (M2-M| : )G 2 (fc 2 ,M2,M s 2 J, (A.5) 
where we assume M 2 = 2B(2m s + m)/3 and where 



(A.6) 



(fc 2 +M 1 2 )(fc 2 +Mf) 
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A. 2 Propagators in SU(4\2) chiral perturbation theory 

We consider the case of the effective theory with Nf = 2 dynamical flavours with two 
degenerate light quarks of mass m and two quenched quarks s and r with the mass of 
the strange quark m s = m Sg = m r = m rg . The dynamical degrees of freedom in the 
effective theory have masses 



Mi 



2Bm, M^ = B(m + m s ), 



Mi 



2Bm, . 



(A.7) 



The propagator G ab (k 2 ) is defined as follows: The Lagrangian in Eq. (|4.3p restricted to 



lowest order and with C\ 



(2) 



£ (2) 



2i> 



- can be rewritten as 



M 



K) 



(A. 



where 



M: 



Ml 

\{Ml + M 2 K ) for 
Ml 



for a = 1, 2, 3 , 



a = 4,5, 6, 7, 9, 10, 11, 12, 16, 17, 18, 19, 25, 26, 27, 28 , 
for a = 8, 13, 14, 15, 20, 21, 22, 23, 24, 29, 30, 31, 32, 33, 34, 35 , 

(A.9) 

and where for our choice of the S'L r (4|2)-generators 



h ah = (V2g a8 + g a15 - g a2A - V2 g a35 )(V2 g b8 + g bl5 - g bM - V2g b35 ) 
l°*> = -(V2g a8 + g al5 + g a2A + V2 g a35 )(V2g b8 + g bl5 + g b24 + V2g b35 ) 

In terms of these matrices the propagator is given by 

G ab (k 2 ) = g ab G 1 (k 2 ,M 2 ) + y ab (M 2 -M 2 K )G 2 (k 2 ,M 2 s ,Ml), 
where G\ and G2 are defined as in Eq. (|A.6|) . 



(A.10) 



B The integral B 21 (q 2 ,M 2 



Here we just quote the integral itself. Its derivation from the sum of the unitary con- 
tribution in Eq. (|4.12[) and the tadpole in Eq. (|4,13p can be found for example in the 
appendix of Ref. [26] : 



B 2 i(q 2 ,M% 
where for q 2 < 4M 2 



i 

48^ 



(B.l) 



B{q 2 ,M 2 ) 



AM 2 



Wit 2 



-2 + v^log ^§±1 



2 ■ 4 
1 „ 21 _i_ I 1 1 
96^ T 960?r 2 IF ~ • • ■ 



(B.2) 
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In the body of this paper we found it convenient to absorb a logarithmic term into the 
expression for B21' 

J B 21 ( / , 2 ,p 2 ,M 2 ) = J B 21 (p 2 ,M 2 )-i i ^log(^) . (B.3) 
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